Crossover Behavior from Decoupled Criticality 



O 

(N 

00 



CO 
-4-3 



i 

o 
o 



> 

On 

m 
cn 

o 
o 



Y. Kamiya, 1 N. Kawashima, 1 and C. D. Batista 2 

1 Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 227-8581, Japan 
2 Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

(Dated: September 29, 2010) 

We study the thermodynamic phase transition of a spin Hamiltonian comprising two 3D magnetic 
sublattices. Each sublattice contains XY spins coupled by the usual bilinear exchange, while spins in 
different sublattices only interact via biquadratic exchange. This Hamiltonian is an effective model 
for XY magnets on certain frustrated lattices such as body centered tetragonal. By performing a 
cluster Monte Carlo simulation, we investigate the crossover from the 3D-XY fixed point (decoupled 
sublattices) and find a systematic flow toward a first-order transition without a separatrix or a new 
fixed point. This strongly suggests that the correct asymptotic behavior is a first-order transition. 

PACS numbers: 75.30.Kz, 05.70.Jk, 05.50.+q, 75.40.Mg 



I. INTRODUCTION 

Geometric frustration can play a decisive role in the 
behavior of magnetic systems. The combination of frus- 
trated geometries with strong quantum fluctuations can 
lead to new quantum states of matter— - — . It has been 
shown recently that novel charge effects in Mott insula- 
tors, such as spin-driven electronic charge density waves, 
or orbital currents, only take place in geometrically frus- 
trated lattices^. Geometric frustration can also reduce 
the effective dimensionality of certain quantum critical 
points^. Finally, it has been known for years that the 
presence of geometric frustration can change the nature 
of certain thermodynamic phase transitions. However, 
it has been also recognized that the nature of the new 
transition can be very elusive for the standard rcnormal- 
ization group treatments^ and may require very sophis- 
ticated numerical approaches^. 

Several quantum magnets comprise two sublattices 
of magnetic ions coupled by a geometrically frustrated 
exchange^. This is for instance the case of a Heisenberg 
antiferromagnet on a body centered tetragonal (BCT) 
lattice^, or a square lattice with nearest- and next- 
nearest-neighbor exchange interaction ;) 10 ' 11 . We are in- 
terested in the regime of inter-sublattice coupling smaller 
than the intra-sublatticc exchange. We will also assume 
that there is a uniaxial easy-plane anisotropy that re- 
duces the Hamiltonian symmetry from 0(3) to 0(2). 
The frustrated nature of the inter-sublattice exchange 
precludes a bilinear coupling between the order parame- 
ters of the two sublattices. The Hamiltonian symmetry 
only allows for an effective biquadratic coupling. Conse- 
quently, if niA and mj are the XY magnetizations at 
wave-vector fco = (tt, -k, 0) of the sublattices (which in 
the BCT lattice case are the even- and odd-numbered 
layers^), the parallel or anti- parallel orientations of 
and m.B correspond to different ground states. The 
Z2 symmetry is broken by selecting one of these two 
statea2»i2r— . The 0(2) XZ2 symmetry breaking also ap- 
pears in the XY model on a triangular lattice^. In this 
case the Z2 broken symmetry corresponds to the two pos- 



sible vector chiral orderings. 

We want to explore the nature of the thermodynamic 
phase transition associated with the 0(2) XZ2 symmetry 
breaking that takes place in several frustrated magnets. 
For this purpose, we will consider classical magnetic mo- 
ments because the quantum character of the spins does 
not affect the nature of the thermodynamic transition. In 
Ref. ©, we used two different approaches to understand 
the effect of the additional Z2 symmetry breaking and 
compared their results. The first approach was a Monte 
Carlo (MC) simulation of the classical spin model on the 
BCT lattice. The second approach was a scaling anal- 
ysis of the Landau-Ginzburg- Wilson (LGW) model that 
preserves the symmetries of the lattice Hamiltonian. A 
single transition with exponents close to those of the 3D 
XY model was obtained from a finite-size scaling (FSS) 
analysis of the MC data£. On the other hand, the scaling 
analysis of the LGW model, 

H LGW = [ d d x\ J2 (\\^^a\ 2 +t\cp a \ 2 +u\<j, a A 
J L a= A,B ^ J 
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indicated that A is a relevant perturbation for the 3D 
XY decoupled fixed point (DFP) located on the u axis 
(u ± 0, A = g = 0)4. Here, </> a = {<f>t,f a ) = A, B) is a 
two-component field representing antiferromagnetic mo- 
ments in even- (a = A) or odd- (a = B) numbered layers, 
and A is the biquadratic coupling between them. These 
results look contradicting at a first glance: although the 
numerical observations can be explained in a consistent 
way by the DFP, this fixed point is nevertheless unstable 
along the A-direction. More specifically, near the DFP, A 
transforms as A' = b Vx X, where y\ = 0.526(8) and b is a 
rescaling factor—. 

The scaling argument implies that there will be a 
crossover behavior from the DFP, provided |A| L Vx > 1 
with L being the system-size^. However, |A| can be quite 
small for the original frustrated spin system because it 
is an effective interaction that arises from second-order 
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perturbation with respect to the ratio between the inter- 
and the intra-layer bilinear exchange couplings^. In ad- 
dition, we could not obtain data for sufficiently large L in 
our previous calculation in Ref. |9| because we simulated 
the original Hamiltonian on the frustrated lattice. Thus, 
the nature of the crossover was left as an open problem. 



II. MODEL AND METHOD 
A. Model 

In this paper, we explore the expected crossover by 
studying an XY spin model on a cubic lattice that is more 
directly related to the LGW effective model than to the 
original Hamiltonian on the BCT lattice. The relevant 
coupling A is explicitly taken into account by considering 
the Hamiltonian model: 



H 
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(i,j),a=A,B 



with J > 0. S ay i (a = A, B) is a classical XY spin at site i 
on the cubic lattice and (i, j) is a pair of nearest-neighbor 
sites. The coefficient A characterizes the amplitude of the 
biquadratic coupling that is expected to drive the system 
away from the DFP. We consider the case A < 0, which 
is experimentally relevant^. No term corresponding to 
the g-term in 7?lgw is explicitly included in H, because 
it is automatically generated when the short wavelength 
modes are integrated out (renormalization process). 

In the ground state, both A and B spins are ferro- 
magnetically ordered and the 0(2) symmetry is bro- 
ken. In addition, their relative phase is locked so that 
Sa,i ■ Sb ;1 = ±1, which causes Z2 symmetry breaking. 
The order parameters associated with these two kinds of 
symmetry breaking arc m = with = L~ d S 0l , 
and (7 = L~ d £\ a\ with tji = Sa,{ ■ Ss,i, respectively. We 
introduce the correlation functions G™ = {Saj ■ Saj) 
and Gfj = (cn&j). (For the definition of m and G™, we 
can use either A or B spins without loss of generality.) 

For very small |A| (A = —0.05, L < 64), we observe 
an apparently continuous transition with exponents of 
the DFP, which is naturally interpreted as the same 
behavior as in the previous MC simulation in Ref. [i| 
However, a more careful FSS analysis reveals the ex- 
pected crossover. We present a numerically obtained 
rcnormalization-group flow diagram of several scaling pa- 
rameters that should be scale-invariant at the second- 
order transitions^. We find that the flow evolves sys- 
tematically from the DFP without a sign of a stable fixed 
point or a separatrix, toward the region where the transi- 
tion is discontinuous. Based on this observation and the 
lack of a stable fixed point in the e-expansion (e = 4 — d) 
around the DFP—, we propose that the correct asymp- 
totic behavior is a first-order transition for any (negative) 
finite value of A. 



B. Method 

The absence of explicit frustration is the main com- 
putational advantage of H relative to the original model 
studied in Ref. [^. This enables us to develop an effi- 
cient cluster MC algorithm based on a minor modifi- 
cation of the embedding method proposed by Wolffii. 
In every update cycle, we choose a unit vector n at 
random. The vector n defines the Z2 transformations 
S A = S A - 2 (S A ■ n) n and S B = -S B + 2 (S B • n) n. 
(The difference by a factor of —1 serves to enhance 
the relaxation of the a modes as compared to apply- 
ing the same mirror-image transformation to the A and 
B spins.) Then, we choose a spin S aj! ; and identify a 
cluster G = {S a; i, Sbj, S Ci fe, . . . } that can be reached 
from S a ^ via probabilistically activated links. The prob- 
ability to activate a link depends on the interaction on 
the link: Pi(S, S') = 1 - min {l, exp [j3J (S - S) • S'] } 
for links with the bilinear exchange and P2(S,S') = 
l-min(l,cxp{|A|/3J (S • S') 2 - (S • S') 2 }) for links 

with the biquadratic coupling. After a cluster is identi- 
fied, we flip it, namely apply the Z2 transformation on 
every spin included in G. It can be easily checked that 
the algorithm satisfies both the detailed-balance and er- 
godicity conditions. 



III. RESULTS 

A. Conventional scaling analysis 

We first present the results for very small |A| with 
|A| L Vx < 1, where we observe an apparently contin- 
uous transition controlled by the DFP. This is natu- 
rally expected from the scaling argument given above 
and basically the same behavior that was observed in 
the frustrated model previously studied in Ref. 0. In 
Fig. DTa), we present the FSS plots of G™ and Gfj at 
the largest distance in a given system where r^j = 
Tij, v = r ijtZ = L/2 (A = -0.05, L < 64). These 
plots are based on the following FSS forms at the DFP—: 
G™{T,L,r. l3 ) ~ i-("+ 1 )/ m ( J L 1 /-( T _ T c ), nj /L) and 
Gfj {T, L, ry) ~ L- 2 ("+ 1 )/ <T ( J L 1 /-( T _ T e ),r«/L) with 
r] = 0.0380(4) and v = 0.68155(27) being the critical 
exponents of the 3D XY models. Using the exponent of 
the DFP, we can also produce reasonable FSS plots for 
the correlation ratios g m and <7cr^, defined by ratios of 
the corresponding correlation functions at two different 
distances = r^y = r^^ = L/2,L/4 [see Fig. [2(b)]. 



However, since the scaling argument shows that the 
DFP is unstable, we conclude that these FSS plots sim- 
ply describe the "pseudo-scaling" behavior, i.e., as long 
as |A| is finite, significant deviations should eventually ap- 
pear in large enough lattices. In other words, we cannot 
conclude that the transition is of second order because 
weak first-order transitions can become practically in- 
distinguishable from continuous transitions in the usual 
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FIG. 1: (Color online) "Pseudo-scaling" behavior observed 
for A = -0.05 (| A| L Vx ^ 0.45 for L = 64) of (a) correlation 
functions at a distance rtj >x = ry )S = ry jJS = L/2 and (b) 
correlation ratios. Here, r\ and v are critical exponents of the 
3D XY model. T c /J ~ 2.2021 is obtained from the crossings 
of dimensionless scaling parameters. 
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FIG. 2: (Color online) Bimodal energy distribution at T : 
for (a) A = -10 (|A| L vx « 43 for L = 16) and (b) A = -2 
(|A| L HA ~ 18 for L — 64). Most error-bars are smaller than 
the symbol-sizes, (c) Peak-to-peak distance of the distribu- 
tion corresponding to the latent heat for the largest L for each 
A. (d) System-size dependence of the peak-to-peak distance. 



FSS analysis for small L. Indeed, for relatively large |A|, 
we find obvious deviations from the DFP. As shown in 
Figs.[2Ja) and (b), the energy distributions near the tran- 
sition show a bimodal structure with increasing depth for 
larger system sizes. This is clear evidence for a first-order 
transition. The peak-to-peak distance gives an estimate 
of the latent heat AE(X). As expected, the first-order 
nature becomes weaker for smaller |A| [see Fig. H[c)]. 



B. Monte Carlo renormalization group analysis 

Given our results for small and large values of | A |, it is 
natural to ask if there is a multicritical point where the 
first-order transition line terminates. The dependence of 
AE (A) on small values of |A| does not provide an efficient 
way of answering this question because larger lattices are 
required to detect smaller values of AE. In what follows, 
we explain our method to investigate the correct asymp- 
totic behavior for very small |A|. Our approach is a sort 
of MC renormalization group analysis^. A similar tech- 
nique was applied, for instance, to the random-bond Ising 
model by Hukushima and it was found that the method 
is very useful to obtain qualitative structure of the phase 
diagram^. 

We consider several dimensionless scaling parameters 
R (A) (such as g m and g a defined above) and introduce 
their L-dependent estimators R (A, L) as the crossings of 
temperature-dependent curves of the parameters for two 
successive system-sizes L and 2L. Because the L — > oo 
limit, R (A), is expected to be scale- invariant and univer- 
sal for a second-order transition, R (A, L) must converge 
to such a universal value if the transition is continuous. 
Consequently, if a multicritical point exists, the "flow" 
structure of R (A, L) should have a separatrix and a stable 
fixed point. Here, the term "flow" refers to the evolution 
of R (A, L) with increasing L. 

In addition to g m and g a , we use as R(X,L) the 
Binder parameters defined by U m = (|m| 4 )/(|m| 2 ) 2 and 
U a = (a 4 ) / (a 2 ) 2 , and the second-moment correlation- 
lengths 2 ^ divided by the system-size £ rn /L and ^ a /L. 
Hence, the entire parameter space is six-dimensional in 
our treatment. The obtained flow diagrams are shown 
in Fig. [3J As can be seen in Figs. [3Jb— d) , we find that 
in the 4D subspace {gm, £,m/L, g ai £,a/L) trajectories of 
the projected flows collapse on an approximately single, 
monotonous curve. Therefore, it turns out to be suffi- 
cient to treat the projected flow in the subspace spanned 
by one of the above four (we choose £ m /L) and the other 
two parameters not included here, namely U m and U a . 

The flow projected onto this (U m ,^ m /L,U a ) subspace 
is shown in Fig. [SJa). The DFP is associated with the 
flows for A = —0.005 or —0.05, because, as is implied by 
the data-collapse in the FSS plots shown in Fig. [TJ with 
such small |A| the effect of the biquadratic perturbation 
is still negligible in the length-scale under consideration. 
The known estimates for the 3D XY universality class are 
U m = 1.2430(5) and £, a /L = 0.5925(2)^. We show the 
point corresponding to these values on the (U m ,£ m / L) 
plane in Fig.[3ta). (Estimates for the other less common 
parameters are not available in the literature as far as we 
know.) The above observation is in good agreement with 
these estimates. 

For larger values of |A| with |A| L yx > 1, the flow clearly 
deviates from the trajectory dominated by the DFP. This 
is a clear sign of the expected crossover. The crossover is 
already evident for A = -0.4 (|A| L» x w 3.1 for L = 48). 
As |A| increases, the flow keeps evolving away from the 
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FIG. 3: (Color online) (a) The flow diagram projected onto 
the (Um,£m/L,U a ) space. Two-dimensional projections are 
also shown. Error bars for U m are shown on the bottom 
plane and those for the other parameters are smaller than the 
symbol size (not shown). System-sizes corresponding to the 
data points in each flow are L = 8, 12, 16, 24 (not for A = 
-0.005) and 32 (only for A = -0.05, -1.4), in order specified 
by arrows attached to the flow lines. The DFP projected on 
the {Um,£m/L) plane is denoted by "XY." U a « 2.40(1) at 
the DFP is estimated by extrapolating the A = —0.005 flow 
and it is shown by a small filled circle on the U a axis, (b)-(d) 
The same flow diagrams projected on the other subspaces. 
The arrows show the overall direction of the flows. 



is relevant for several frustrated magnets near their ther- 
modynamic phase transitions. Such crossover results in 
a weakly first-order phase transition. Our calculation 
also shows that it will be very difficult to observe such a 
first order transition with standard experimental meth- 
ods as long as the frustrated inter-layer coupling is small 
in comparison with the intra-layer exchange. This is in- 
deed the case of BaCuSi2 0g2i as discussed in Ref. 0. In 
other words, although the correct asymptotic behavior 
is the first-order transition, the thermodynamic behavior 
will be dominated by the 3D XY DFP in a broad re- 
gion near the transition. The true discontinuous nature 
of the transition can be observed in a very narrow region 
near the transition point that could easily be beyond the 
experimental precision in most cases. Nevertheless, the 
first order transition should be observable for frustrated 
magnets with |A| of order one. In such cases, the 3D XY- 
like behavior beyond a certain distance from the transi- 
tion point will be finally interrupted by the fluctuation- 
induced first-order transition. 

A value of |A| of order one is indeed realized in the frus- 
trated spin model that has been proposed for describing 
the iron based superconductors LaFeAs(Oi- a; F a ;) 22 ' 23 . 
According to our result, such a model should exhibit 
a single weakly first-order transition to the broken 
0(2) XZ2 phase in presence of a strong magnetic field 
(the field is required to induce effective 0(2) magnetic 
moments). The stacked triangular antiferromagnetic 
compounds^ are other physical realizations of the effec- 
tive model considered here [Eq.([2])]. Similarly, we predict 
a single weakly first-order phase transition to take place 
in these systems in presence of a strong magnetic field, 
which is in agreement with recent investigation s 8 ' 25 . 



DFP without a stable fixed point or a separatrix. Note 
that we have already shown clear evidence of a first-order 
transition for A = — 2 [Fig. Hfb)]. This indicates that 
the observed crossover eventually leads to the first-order 
transition. 

While the numerical evidence in finite systems is al- 
ways insufficient for very small |A|, we take the numer- 
ical result presented above as a strong evidence for the 
first-order character of arbitrary small |A|. This conclu- 
sion is also supported by the epsilon expansion analysis 
of -Hlgw around the DFP—: the result obtained by ex- 
panding the Hamiltonian to O(e) is most naturally ex- 
plained as the lack of a separatrix fixed point, suggesting 
a fluctuation-induced first-order transition. 



IV. SUMMARY 

To summarize, we have established the crossover be- 
havior from the 3D XY DFP for an effective model that 
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